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ABSTRACT 

We use deep far-IR, submm, radio and X-ray imaging and mid-IR spectroscopy to explore the nature 
of a sample of Spiteer-selected dust-obscured galaxies (DOGs) in GOODS-N. A sample of 79 galaxies 
satisfy the criteria R — [24] > 14 (Vega) down to S24 > 100 /xJy (median flux density S24 = 180 /xJy). 
Twelve of these galaxies have IRS spectra available which we use to measure redshifts and classify 
these objects as being dominated by star formation or active galactic nuclei (AGN) activity in the 
mid-IR. The IRS spectra and Spitzer photometric redshifts confirm that the DOGs lie in a tight 
rcdshift distribution around z ~ 2. Based on mid-IR colors, 80% of DOGs are likely dominated 
by star formation; the stacked X-ray emission from this sub-sample of DOGs is also consistent with 
star formation. Since only a small number of DOGs are individually detected at far-IR and submm 
wavelengths, we use a stacking analysis to determine the average flux from these objects and plot a 
composite IR (8-1000 /zm) spectral energy distribution (SED). The average luminosity of these star 
forming DOGs is Lir ~ 1 x 10 12 Lq. We compare the average star forming DOG to the average bright 
(5*850 > 5mJy) submillimeter galaxy (SMG); the S24 > 100 ^iJy DOGs are 3 times more numerous 
but 8 times less luminous in the IR. The far-IR SED shape of DOGs is similar to that of SMGs 
(average dust temperature of around 30 K) but DOGs have a higher mid-IR to far-IR flux ratio. The 
average star formation-dominated DOG has a star formation rate of 200M Q yr _1 which, given their 
space density, amounts to a contribution of 0.01 M Q yr -1 Mpc~ 3 (or 5-10%) to the star formation 
rate density at z ~ 2. We use the composite SED to predict the average flux of DOGs in future 
Herschel/PACS 100 /im and SCUBA-2 450 surveys and show that the majority of them will be 
detected. 

Subject headings: galaxies: active — galaxies: evolution — galaxies: starburst — infrared: galaxies 
— submillimeter — X-rays: galaxies 



1. INTRODUCTION 

Large extragalactic surveys with the Spitzer Space 
Telescope (Werner et al. 2004) have revealed many 
high redshift objects which are bright in the mid- 
infrared (mid-IR) and have red mid-IR to optical colors 
(e.g. Houck et al. 2005; Yan et al. 2007). The selection in 
the mid-IR indicates that these objects are dusty; how- 
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ever, without mid-IR spectra it is not clear whether the 
dust is heated by active galactic nuclei (AGN) or star 
formation activity or both. Furthermore, extrapolating 
from mid-IR to total IR luminosity is uncertain without 
good constraints spanning the far-IR dust peak (e.g. Pa- 
povich et al. 2007; Daddi et al. 2007). 

A sample of IR-luminous Dust Obscured Galaxies 
(DOGs), selected to have very red i?-[24] color, was re- 
cently presented by Dey et al. (2008, hereafter D08). 
From spectroscopic observations DOGs were found to 
have a tight redshift distribution around z ~ 2, very sim- 
ilar to that of the submillimeter selected galaxies (SMGs, 
e.g. Chapman et al. 2005). The space density and clus- 
tering of DOGs are also comparable to those of the bright 
(5850 > 6 mJy) SMGs, suggesting that these two popu- 
lations might be associated (e.g., in an evolutionary se- 
quence; D08, Brodwin et al. 2008). 

The origin of the bolometric luminosity in DOGs is 
uncertain. Using a similar i?-[24] selection and an addi- 
tional R-K color criterion, Fiore et al. (2008, hereafter 
F08) selected a sample of faint DOGs (5*24 > 40 /iJy) in 
the CDF-S and concluded, based on their stacked X-ray 
spectrum, that 80% are Compton-thick AGN. 

In order to test the relationship between SMGs, DOGs, 
and the F08 Compton-thick AGN, one needs multi- 
wavelength data including near-IR, far-IR and submil- 
limeter observations of a large sample of DOGs. In 
this paper we use the deep multi-wavelength data in the 
GOODS-N field (Giavalisco et al. 2004) to study a sam- 
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pie of faint DOGs in order to put constraints on their in- 
frared luminosities, determine the role of AGN and star 
formation activity in these systems and compare them 
with SMGs. Our goal is to improve our understanding 
of the role of DOGs in massive galaxy evolution. 

All magnitudes in this paper use the Vega system un- 
less otherwise noted. We assume a standard cosmology 
with H = 72kms- 1 Mpc~\ fi M = 0.3 and Q A = 0.7. 

2. DATA AND SAMPLE SELECTION 

Using the Subaru i?-band (Capak et al. 2004) and 
Spitzer 24 /zm (Chary et al. in preparation) data in 
GOODS-N we select 79 galaxies which satisfy the DOGs 
criteria of i?-[24]> 14 (i.e. S 2 a/S r > 1000) and S* 24 > 
100 /iJy. The 24 /im fluxes are measured by fitting the 
point spread function using the IRAC positions as pri- 
ors (see Chary et al. in preparation for details), and the 
i?-band photometry is measured in 3 arcsecond diameter 
apertures since this gave optimal SNRs. This GOODS- 
N sample is significantly deeper than the Bootes DOGs 
sample of D08, which was limited by S24 > 300 /xJy; 
there are only 13 DOGs in GOODS-N above this flux 
cut. The median 24 /im flux density for the 79 DOGs in 
GOODS-N is 180 /iJy. 

To understand the nature of DOGs we exploit the deep 
multi-wavelength data available in GOODS-N. We mea- 
sure the if -band magnitude of each DOG in the new 
WIRCAM/CFHT if-band images (Lin et al. in prepara- 
tion), and the B and z magnitudes from the Capak et 
al. (2004) Subaru images in matched circular apertures. 
For comparison we also measure the R and K magni- 
tudes for the sample of SMGs in GOODS-N (Pope et 
al. 2006). The IRAC fluxes of all DOGs are available 
from the deep Spitzer legacy images (Dickinson et al. in 
preparation), where we use 4 arcsecond diameter aper- 
tures (with the appropriate aperture corrections applied) 
for photometry. 

Many DOGs are not individually detected in the X- 
ray, far-IR, (sub)mm, and radio (see Table [T]) and so we 
must rely on stacking analyses. We use the Chandra X- 
ray Observatory 2 Msec image (Alexander et al. 2003) 
to identify DOGs which are formally detected and also 
to perform a stacking analysis of the undetected sources. 
MIPS 70 and 160 /im data were reduced following the 
techniques of Frayer et al. (2006a,b). The MIPS 70 and 
160 /im images are searched for detections as well as used 
for stacking analysis (see Huynh et al. 2007 for details). 
In the (sub)mm, we use the 850 /im Submillimetre Com- 
mon User Bolometer Array (SCUBA) super-map (Bo- 
rys et al. 2003; Pope et al. 2005) and the 1.2 mm Max- 
Planck Millimetre Bolometer (MAMBO) map (Greve et 
al. 2008). For (sub)mm stacking we perform a variance- 
weighted average due to the variable noise levels across 
the maps and stack on an image which has the point 
sources removed so as not to bias the result. At the 
longest wavelengths, we use the 1.4 GHz VLA A+B-array 
map and catalogs of Morrison et al. (in preparation). 
Stacking in the radio and far-IR is performed by stack- 
ing images centered on each source and therefore we can 
check to make sure that the resulting point spread func- 
tion is as expected. For the stacking analyses in the 
X-ray, MIPS, SCUBA, MAMBO and radio images we 
perform Monte Carlo simulations at random positions to 
determine the error and significance of the stacked result. 
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Fig. 1. — S24/SR as a function of R— K color for DOGs and 
SMGs in GOODS-N. The dashed line indicates the DOGs selection 
criteria while the dotted line shows the additional color constraint 
in Fiore et al. (2008). At least 90% of DOGs and 20% of SMGs 
meet the Fiore et al. (2008) criteria (30% of SMGs are classified as 
DOGs). 

Several Spitzer IRS spectroscopy programs have tar- 
geted GOODS-N; we find twelve DOGs have spectra 
available from GO-20456 (10 sources, Pope et al. 2008; 
Murphy et al. in preparation) and GO-20733 (2 sources, 
J. Van Duyne private communication). Details of the 
IRS data reduction can be found in Pope et al. (2008). 

3. RESULTS 
3.1. Comparison with other selection criteria 

Fig.[T]shows the R-K-S24 color-color plot (cf., F08) for 
the DOGs and the GOODS-N SMGs (Pope et al. 2006). 
While SMGs are not constrained to any region of the 
diagram, the DOGs almost all satisfy the additional R — 
K criterion from F08 (dotted line; F08 conclude that 
80% of galaxies in this region are Compton-thick AGN). 
This is somewhat unexpected, since Fig. 3 of F08 shows 
just as many S24/SR > 1000 galaxies on either side of 
R — K = 4.5. However F08 push even deeper (S24 > 
40/iJy) and it is clear that the fainter 24 /im samples 
contain a higher fraction of bluer (in R — K) objects 
(F. Fiore, private communication). Fig.[T]also shows that 
30% (> 20%) of SMGs meet the DOGs (F08) criteria, 
respectively. It is difficult to assess the fraction of DOGs 
which arc submm detected, due to the highly varying 
noise levels in the SCUBA map of GOODS-N (Borys 
et al. 2003), but looking only at the low noise regions 
(< 2.3 mJy RMS) of the map it appears that around 
30% (7/24) of DOGs are formally detected (> 3.5cr) at 
850 /im. In contrast, only 2/73 DOGs are coincident with 
detections in the GOODS-N MAMBO map, down to a 
1.2mm depth of - 0.8 mJy RMS (Greve et al. 2008). 
The difference between the submm and mm detection 
rates of DOGs could be because the millimeter selection 
picks out either higher redshift or cooler objects than the 
submillimeter (see discussion in Greve et al. 2008). 

Another widely-used selection criteria for high redshift 
galaxies is the BzK color-color plot (Daddi et al. 2004). 
Fig. E] shows the BzK plot for the DOGs in GOODS- 
N. All DOGs with a > 3er detection in K are plotted 
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Fig. 2.— BzK color-color plot for DOGs in GOODS-N. The 
diagonal line separates active galaxies at z = 1.4—2.5 (above line, 
Daddi et al. 2004). Squares and triangle are the GOODS-N DOGs 
where the different symb ols show the mid-IR SF and AGN classified 
DOGs (see Section l3,3j l, 

(75/79); many of these are only limits at B and z. While 
almost all DOGs are consistent with the BzK selection, 
many would not be included in robust active BzK sam- 
ples since they are not detected in B and/or z. Roughly 
12% of BzK galaxies down to S24 > 100 /iJy will satisfy 
the DOGs criteria. The overlap between DOGs and BzK 
galaxies is an interesting topic, and will be explored in 
more detail in a future paper (Meger et al. in prepara- 
tion). The position of the DOGs within the BzK plot 
is our first clue that these are star-forming galaxies at 
z = 1.4-2.5 (Daddi et al. 2004); the redshifts, star for- 
mation and AGN activity of DOGs will be discussed in 
the following sections. 

3.2. Redshifts 

A redshift distribution peaking around 2 is expected 
for the DOGs selection criteria since the 24 /jm detec- 
tion will limit the sample to sources below z ~ 3 and 
the red mid-IR to optical color allows only for sources 
which are very faint in the optical, which weeds out many 
z < 1 sources. D08 presented a redshift distribution for 
bright (S24 > 0.3 mJy) DOGs based on spectroscopically 
measured redshifts, resulting in a Gaussian distribution 
with z=2.0 and <r(z)—0.5. D08 also note no obvious 
dependence of redshift on 24 /jm flux within their sam- 
ple. F08 presented a photometric redshift distribution 
for galaxies which meet their two color selections down 
to S24 > 40 /iJy, and it is remarkably similar to that of 
the D08 DOGs. This suggests that our S 2 i > 100 /xJy 
sample of DOGs will also have a similar redshift distri- 
bution to that of the brighter DOGs in D08. 

The redshifts for the 12 GOODS-N DOGs with IRS 
spectra range from 1.6-2.6 with a peak at 2, consistent 
with the spectroscopic redshift distribution from D08. 
However, 10/12 of the GOODS-N DOGs with IRS spec- 
tra are brighter than 5*24 > 300 /iJy. In addition to the 
12 DOGs with IRS spectra, only 1 additional DOG has 
an optical spectroscopic redshift from the ~ 3000 spec- 
troscopic redshifts available in GOODS-N emphasizing 
how faint these galaxies are in the optical. This DOG is 



Fig. 3.— Redshift distribution for GOODS-N DOGs using the 
Pope et al. (2006) photometric redshift estimator (thin solid dis- 
tribution). The thick solid histogram shows the DOGs with IRS 
spectroscopic redshifts. The dotted curve is the Gaussian fit to the 
photometric redshifts which gives < z >= 2.0 ± 0.3. 



the unusual Waddington et al. (1999) dusty radio galaxy 
at z — 4.424; the highest confirmed redshift 24 /im source 
in GOODS-N. 

Although most DOGs are too faint in the optical to 
yield accurate traditional optical photometric redshifts, 
we can use the IRAC photometry as a rough estimate 
of the redshift since at z = 2 these channels sample the 
1.6 fim stellar bump (e.g. Simpson & Eisenhardt 1999, 
Sawicki 2002). Recently, Farrah et al. (2008) presented 
a study of l bump2' sources defined as 53.6 < S4.5 > S5.8 
and S4.5 > 58.o and found them to lie in a tight redshift 
distribution, < z >— 1.71±0.15. Similarly, we can select 
'bump3' sources (54.5 < S5.8 > Ss.o and S3. 6 < S5.8) 
which should lie around z ~ 2.5. Of the 79 DOGs 
in GOODS-N, 62 (78%) satisfy either the l bump2' or 
'bump3' criteria which places them between z — 1.3-2.9; 
this is consistent with the spectroscopic redshift distribu- 
tion of DOGs in D08. The presence of a stellar bump in 
most DOGs also indicates that they are not dominated 
by AGN emission in the rest-frame near-IR; this will be 
discussed further in the next section. 

As an additional test of the redshifts of DOGs, we used 
the independent photometric-redshift estimate derived in 
Pope et al. (2006). Equation 2 of Pope et al. (2006) pro- 
vides a simple empirical relation between the redshift and 
the IRAC and MIPS 24 photometry which has been 
tuned to SMGs with spectroscopic redshifts. The appli- 
cation of this relation to the DOGs may be justified on 
the basis that both DOGs and SMGs are dust obscured 
populations of galaxies, and may therefore have similar 
SEDs. Indeed, the median IRAC flux densities of the 
DOGs and SMGs in GOODS-N are almost identical. 

We first tested this photometric redshift technique us- 
ing the Bootes DOGs with spectroscopic redshifts (D08). 
We removed the DOGs which show a power-law in IRAC 
since this method relies on the presence of the stellar 
bump to estimate the redshift. This sub-sample of 29 
Bootes DOGs has a median spectroscopic redshift of 1.9 
(interquartile range 1.6-2.1) and the photometric red- 
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shifts using Equation 2 of Pope et al. (2006) have a distri- 
bution with a median of 1.8 (interquartile range 1.4-2.1). 
Comparing the spectroscopic and photometric redshifts 
for the individual Bootes DOGs we measure a scatter of 
er(Az/(l + z)) — 0.3 with no obvious biases in the pho- 
tometric method. This photometric redshift estimator 
appears to do a reasonable job in estimating the redshift 
distribution of DOGs, although on an individual basis 
the redshifts are still quite uncertain. The redshift dis- 
tribution for the GOODS-N DOGs using Equation 2 of 
Pope et al. (2006) is shown in Fig. [3] (thin solid distribu- 
tion) along with the subset of DOGs with IRS spectro- 
scopic redshifts (thick solid histogram). The photometric 
redshift distribution confirms that most of our GOODS- 
N DOGs sample lies between z = 1.5-2.5. 

For the rest of this paper, we focus on studying the 
multi-wavelength properties assuming the average DOG 
is at z — 2. 
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3.3. AGN and star formation activity 

The IRS spectra also help to quantify the contribu- 
tion from AGN and star formation (SF) activity to the 
mid-IR luminosity of DOGs. We use the same spectral 
decomposition described in Pope et al. (2008) and clas- 
sify an object as AGN-dominated if > 50% of the mid- 
IR luminosity is coming from the continuum component. 
We find that of the 12 DOG IRS spectra, 6 are AGN 
dominated, and 6 are SF dominated. As was done in 
Pope et al. (2008, see also Ivison et al. 2004), we plot 
these DOGs on a Spitzer color-color plane (Fig. [4|) and 
find that the IRS SF and AGN dominated sources (open 
squares and diamonds, respectively) separate very nicely 
at Ss.o/Si.5 — 2. This color cut for separating the SF 
and AGN dominated DOGs is consistent with the sim- 
ulations of dusty z = 2 galaxy templates by Sajina et 
al. (2005), and also with the 'bump' criteria of Farrah et 
al. (2008). 

Based on the IRS spectral results, we use this color 
cut (shown by the dashed line in Fig. [4j to separate the 
GOODS-N DOGs into SF and AGN dominated classes, 
and find that 80% fall into to the SF-dominated class. 
Interestingly this is the same fraction found using the 
same diagram for SMGs (Pope et al. 2008). The fainter 
DOGs (5*24 = 100 — 300 /iJy) contain a higher fraction 
of SF DOGs (60/66 ~ 90%) than the 5 24 > 300/iJy 
DOGs (5/13 - 40%). The median 24 /zm flux densities 
for the SF and AGN DOGs are 175 and 310 ^Jy, respec- 
tively. This is consistent with the previously noted trend 
where AGN contribution increases with 5*24 (e.g. Brand 
et al. 2006; D08). 

We note that a few of the AGN DOGs have very red 
S24/SS.0 colors inconsistent with the redshifted Mrk 231 
spectrum. At z = 2, the 24 fim flux can be enhanced by 
the 7.7 fim PAH feature on top of the AGN continuum 
leading to a higher S24/SS.0 color. Alternatively, the 
AGN DOGs could have higher SW/Ss.o colors because 
they contain a heavily obscured AGN; adding more ob- 
scuration to a power-law AGN will decrease S$.o more 
than S*24 leading to redder colors (Sajina et al. 2007). 
We find that the two AGN DOGs which have both red 
S24/SS.0 colors as well as IRS spectra show a continuum 
dominated mid-IR spectrum and have no X-ray detection 
in the deep Chandra images indicating that they harbor 



Fig. 4. — Spitzer color-color diagram used to separate the SF 
and AGN dominated DOGs. Large circles are the DOGs with 
S24, > 300 ^Jy and small circles are DOGs with S24 = 100—300 /iJy. 
Open squares and diamonds are DOGs with IRS spectra classified 
as SF and AGN dominated, respectively. Based on DOGs with 
IRS spectra we classify SF DOGs as having Sg.o/S^.s < 2 (vertical 
dashed line); 80% of DOGs satisfy this criterion. The colors of the 
M82 (starburst galaxy, Forster Schreiber et al. 2003) and Mrk 231 
(AGN-dominated ULIRG, Rigopoulou et al. 1999) as a function of 
redshift are plotted as the dotted and dashed curves, respectively, 
with the numbers corresponding to the redshift. 

Compton-thick AGN 14 . This suggests that this Spitzer 
color-color plot might also be useful to further separate 
obscured and unobscured AGN via the S24/SS.0 color 
although much larger samples are needed to test this. 

We can also use the deep X-ray imaging to investi- 
gate the presence of X-ray emitting AGN in the DOGs 
sample. We find a higher fraction of X-ray detections 
in the AGN DOGs than in the SF DOGs (Table [J). 
7/7 AGN DOGs with X-ray detections have an effec- 
tive photon index of T < 1.0 and an X-ray luminosity of 
> 10 42 ergs _1 (assuming z ~ 2) which indicates that the 
X-rays are coming from obscured AGN emission (Alexan- 
der et al. 2005, 2008); only 3/7 SF DOGs with X-ray 
detections satisfy these criteria. Stacking the X-ray un- 
detected DOGs in the central 6.5 arcminute (radius) re- 
gion of the X-ray image (see Section 4.2 of Alexander 
et al. 2008 for details) gives a detection in the full and 
soft bands, but not in the hard band (Table [2|). This 
implies a hardness ratio (H/S) of < 0.8 (3cr upper limit) 
which differs from the value of 1.3 found in F08 although 
their definition of the hard and soft band is slightly dif- 
ferent from ours (Table [2]). Converting the F08 values 
to the same bands that we use, their stacking analysis 
gives a hardness ratio of 0.8 (r = 1.0). Our 3cr upper 
limit on the hardness ratio of SF DOGs is within the 
error bar of the measured value of F08. A number of 
factors could lead to a difference in the X-ray stacking 
results including the area of the X-ray image used in 
the stacking analysis, the sample size (F08 are stacking 
111 objects), and of course the limiting depth at 24 /im. 
Stacking the mid-IR classified SF and AGN DOGs sep- 
arately we obtain similar results and conclude that the 

14 One of these sources is also an SMG (source CI, Pope et 
al. 2008; see also Alexander et al. 2008). 



The nature of DOGs 



5 



TABLE 1 10"' 

Multi-wavelength detections of DOGs 



Type 


N 


1.4GHz a 


70Atm b 0.5- 


SkeV c 


SF DOGs d 


65 


43 





7 


AGN DOGs 


12 


11 


5 


7 


No IRAC e 


2 











All DOGs 


79 


54 


5 


14 



a ct ~ 5.3 /^Jy (Morrison et al. in preparation) 
b cr ~ 0.6 mjy (Huynh et al. 2007) c a ~ 2.4 X 
10- 17 ergcm- 2 s- 1 (Table 9 of Alexander et al. 2003) 
d Sg/Si.s < 2 c These DOGs are outside the uniform 
IRAC coverage. 



stacked X-ray hardness ratio of the two sub-samples can- 
not be distinguished within our uncertainty. For the SF 
DOGs we estimate L .5-8kcV = 5.7 x 10 41 ergs _1 and 
£2-iokcV = 3.9 x 10 41 ergs _1 from the 1.5-6 keV ob- 
served luminosity assuming T = 1.8. This X-ray lumi- 
nosity is several orders of magnitude lower than that of 
the bright SMGs (Alexander et al. 2005). Depending 
on how much of the X-ray emission is coming from high 
mass X-ray binaries (HMXBs), we convert the X-ray lu- 
minosity of DOGs to SFR and obtain 82-390 Mgyr" 1 
(Bauer et al. 2002 and Persic et al. 2004, respectively). 
With the small sample of DOGs in GOODS-N, the hard 
X-rays are poorly constrained and we cannot rule out an 
X-ray emitting AGN in some DOGs. The high X-ray lu- 
minosities of the AGN DOGs are characteristic of AGN 
activity (e.g. Alexander et al. 2005). However, contrary 
to what is implied by F08, most SF DOGs have X-ray 
emission which is consistent with what is expected from 
star formation (see Section [3^)1 . 

3.4. IR SED 

In order to constrain the IR luminosities of DOGs we 
need data in the far-IR and submm. We have shown 
that only a handful of DOGs are detected in the far- 
IR and/or submm, so we perform a stacking analysis to 
obtain the average SED shape and Lir for DOGs. For 
the following analysis we focus only on the SF DOGs 
since there are not enough AGN DOGs in GOODS-N to 
obtain a representative composite SED. 

Fig.[5]shows a composite SED for the SF DOGs, where 
we plot the stacked average or median values of S24, <SVo, 
S160, S$5o and S1200 (Table [3]) 15 assuming the average 
DOG is at z = 2. The solid curve is the best-fit SED 
to the solid points where we have fit scaled Chary & El- 
baz (2001, CE01) templates with additional extinction 
from the Draine (2003) models (see Pope et al. 2006 for 
further details on the SED fitting). The dotted curve is 
a composite SED for bright (Ssso > 5mJy) SMGs with 
mid-IR spectra (Pope et al. 2008), scaled down by a fac- 
tor of 8, which matches the DOGs SED in the far-IR and 
submm, but is too faint at 24 fim and 70 /xm observed. 
This is consistent with the results of Sajina et al. (2008) 
where bright (5*24 ^ 1 nrJy) high redshift ULIRGs are 
rarely detected at (sub)mm wavelengths (see also Lutz 
et al. 2005). The excess emission in the mid-IR rela- 
tive to SMGs in the sub-sample of strong PAH sources 
from Sajina et al. (2008) accounts for 30% of the total 

15 Note that not all DOGs are used in the stack at each wave- 
length, since the samples for stacking were dependent on the cov- 
erage and depths of the multi- wavelength maps. 




Rest wavelength (|lm) 

Fig. 5. — Composite SED of SF DOGs: we fit the average 
fluxes (solid points) of SF DOGs to the CEOl+Drainc models (solid 
curve). The dotted curve is a normalized (divided by a factor of 8) 
composite SED for SMGs (Pope et al. 2008), and the dashed curve 
is the scaled SMG composite with additional hot (T = 350 K) dust. 
The short horizontal lines indicate the 5<r depths of the planned 
deep surveys at 100 and 450 /an with Herschel/PACS and SCUBA- 
2 and show that the majority of DOGs will be detected by these 
surveys. 

IR luminosity. If we add a hot (T = 350 K) dust compo- 
nent to the SMG composite then we obtain a good fit to 
the SF DOGs (dashed curve). The additional hot dust 
component accounts for less than 10% of the total IR lu- 
minosity and could be due to SF or AGN activity (Tran 
et al. 2001). With IRS spectra existing for only a small 
subset of the brightest sources in our sample we cannot 
say whether this mid-IR excess is due to hot dust or en- 
hanced PAH emission (the SNR of the IRS spectra is not 
high enough to differentiate the PAH equivalent widths 
between the SMGs and the DOGs). However, regardless 
of the source of the excess, our best-fit SED shows that 
the total IR luminosity in these objects is dominated by 
the cold dust component, presumably fueled by star for- 
mation. The 8-1000 /im total Lir for both the solid curve 
and the dashed curve is 1 x 10 12 L© and the average dust 
temperature is ~ 32 K. The IR luminosity implied by the 
median radio flux of DOGs (<~ 20 /iJy) and the radio-IR 
correlation is consistent with this estimate. 

In order to estimate the uncertainty in Lir for the SF 
DOGs we perform 1000 Monte Carlo simulations where 
we randomly sample each data point assuming a Gaus- 
sian with mean equal to the stacked value and a equal 
to the uncertainty in this value. We also include the 
redshift uncertainty assuming a Gaussian redshift distri- 
bution centered on 2 with a = 0.3 (Fig. [3]). The resulting 
distribution of Lir is (1.1 ± 0.5) x 10 12 L , where the er- 
ror is the la uncertainty 16 . This Lir for the SF DOGs 
implies a star formation rate of 200M Q yr _1 , using the 
Kennicutt (1998) relation. This is consistent with the 
range of estimates from the X-ray emission given that X- 
ray SFR relation depends strongly on the relative contri- 
butions from HMXBs and LMXBs to the X-ray emission 
(e.g. Persic et al. 2004). 

1B If we assume a = 0.5 for a Gaussian redshift distribution 
centered on 2 then we get L IR = (1.1 ± 0.7) X 10 12 Lq. 
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TABLE 2 

X-RAY STACKING OF DOGS 



Type 


N 


Counts (10- 6 s" 1 ) a 




H/S 


r 


Flux (10 


- 17 cgs) b 




Luminosity (10 42 ergs 


-1JC 






0.5-8 keV 0.5-2 keV 2- 


-8keV 






0.5-8 keV 0.5- 


2 keV 2- 


-8keV 


1.5-24 keV 1.5-6 keV 6 


-24 ke V 


All DOGs 


30 d 


4.0(5.0<t) 2.2(5.8cr) 


< 1.7 


< 0.8 


> 1.0 


4.7 


1.1 


< 3.9 


1.4 0.33 


< 1.1 


SF DOGs 


28 


3.3(3.9ct) 2.1(5.5<7) 


< 1.8 


< 0.8 


> 1.0 


3.9 


1.0 


< 3.9 


1.2 0.31 


< 1.2 


AGN DOGs 


2 


12.8(3.0(7) < 5.2 


< 8.2 


n/a 


n/a 


15 


< 2.6 


< 19 


4.4 < 0.77 


< 5.6 



a We list the counts in each band if > 3er otherwise we list the 3<r upper limit. Assuming T ■ 
restrict the stacking analysis to the central 6.5 arcmin region of the Chandra images. 



1.4. c Rest-frame X-ray luminosity assuming z — 2. 
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TABLE 3 
Average flux densities of 
GOODS-N SF DOGs 



Wavelength (/im) 


Flux (mjy) 


24 


0.17±0.06 


70 


0.44±0.11 


160 


6.6±2.5 


850 


0.95±0.30 


1200 


0.61±0.10 



The short horizontal lines in Fig. O indicate the 5er lim- 
its of the deepest surveys to be done with Herschel/PACS 
at 100 /mi (Pilbratt 2001) and JCMT/SCUBA-2 at 
450 /Lim (Holland et al. 2006); the majority of DOGs will 
be detected. These surveys will put constraints on the 
infrared luminosities and dust temperatures of individual 
galaxies without the need for stacking. 

4. DISCUSSION 

The average Ltr derived for the DOGs using submm 
and far-IR measurements is a factor of ~ 4 times smaller 
than that calculated in D08 using a conversion from S24 
observed (aka Lg rest = uL^g^m) to Ltr. Part of this 
is because our average L% is lower, since we push three 
times deeper at 24 /im, and part is because of the as- 
sumed conversion between mid-IR and total IR lumi- 
nosity. Based on our best-fit SED, we calculate that 
Lib./ — 7 (quartile range from Monte Carlo simula- 
tions is 5-10) which is within the lower range of con- 
version factors assumed in D08 (L\r/ L$=5-15). On 
the other hand the SMGs from Pope et al. (2008) with 
5*850 > 5mJy have an average Lift/ig of ~ 20. The 
conversion between 24 /im flux and Lir is uncertain for 
high redshift galaxies since it relies on local galaxy tem- 
plates. Observations of high redshift ULIRGs indicate an 
evolution in SED shapes from local ULIRGs (e.g. Pope 
et al. 2006, 2008; Rigby et al. 2008). The overestimate 
of iiR (and SFR) using only S24 has been noted for 
bright high redshift ULIRGs (e.g. Papovich et al. 2007, 
Daddi et al. 2007). It is clear from our full SED fits 
for DOGs and SMGs that a uniform conversion cannot 
be applied to all high redshift ULIRGs and there must 
be additional parameters, other than mid-IR luminosity, 
which are needed to determine the total Ltr. Progress in 
determining these parameters will be facilitated with fu- 
ture wide-field far-IR and submm surveys with Herschel 
Space Observatory and SCUBA-2, for example. 

As discussed in D08, DOGs with S24 > 300 /im have 
similar surface densities to SMGs with Ssso > 6mJy 
(Coppin et al. 2006). The number density of DOGs in 
GOODS-N down to 5 24 = 100 /tJy is 0.5arcmi n - 2 , 6 
times more than the shallower sample in D08. Given the 
average SFR from L IR for the SF DOGs (200M Q yr- 1 ), 



we calculate a SFRD of = 0.01 Moyr^Mpc" 3 at z = 2. 
Depending on the value adopted for the total SFRD at 
z = 2 (e.g. Chary & Elbaz 2001; Caputi et al. 2007), SF 
DOGs contribute 5-10% of the total SFRD at z = 2. 
This result appears to conflict with D08 who suggest 
that (S 24 > 300 /tJy) DOGs contribute 25% of the total 
IR luminosity density at z = 2. However, if we remove 
the AGN DOGs from the D08 sample and use the lower 
conversion factor of Lir/L s ~ 7, this 25% becomes 7% 
which is consistent. The AGN DOGs will further con- 
tribute to the SFRD, however the small number of AGN 
DOGs in GOODS-N does not allow us to put constraints 
on their average Lm and SFR. While the simple DOGs 
selection can be used in many of the deep Spitzer surveys 
to isolate large samples of high redshift ULIRGs, this se- 
lection alone does not resolve the bulk of the SFRD at 
z = 2. For comparison, bright (Ss5o > 5mJy) SMGs con- 
tribute 0.02M Q yr- 1 Mpc" 3 , with fainter (S 850 > 2mJy) 
SMGs contributing 0.05 MQyr^Mpc" 3 at z = 2 (Wall 
et al. 2008). In addition to the DOG samples, several 
different selection criteria using Spitzer data have been 
presented in the literature to isolate ULIRGs at high red- 
shifts (e.g. Yan et al. 2007; Farrah et al. 2008), although 
none of these select a population as numerous as the 
DOGs. Table [1] shows that 70% of DOGs are detected 
in the radio. Since DOGs by definition are faint in the 
optical, these radio-detected DOGs would be similar to 
the optically-faint radio galaxies (OFRGs) discussed in 
Chapman et al. (2002, 2004). 

SMGs are thought to be galaxies in the early stage 
of a massive merger (e.g., Conselice et al. 2003; Pope 
et al. 2008; Tacconi et al. 2008), and D08 propose that 
bright DOGs might be a later stage in the merger. In 
support of this, we find that ~ 30% of SMGs meet the 
DOGs criteria and all 3 SMGs which have > 50% AGN 
contribution in the mid-IR (Pope et al. 2008) are in our 
sample of AGN DOGs. This implies that DOGs selection 
preferentially picks up the more AGN dominated SMGs, 
although these are among the most luminous DOGs. The 
average SF DOG shows additional mid-IR emission com- 
pared to the normalized SMG SED, which may be en- 
hanced PAH emission or hot dust heated by an AGN or 
star formation. Regardless of the source of the mid-IR 
excess emission (which accounts for < 10% of the total 
IR luminosity), the average Ltr and X-ray luminosity 
of the SF DOGs is several times less than that of most 
SMGs indicating that the average DOG is not likely to 
evolve from SMGs. Fig. [2]shows that most DOGs satisfy 
the BzK selection; while they have ULIRG-like luminosi- 
ties BzK galaxies are thought to be forming stars con- 
tinuously over longer timescales and do not necessarily 
require a major merger as as catalyst for star forma- 
tion (Daddi et al. 2008). In summary, we remind the 
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reader that this analysis is focussed on the average prop- 
erties of DOGs in GOODS-N; while the average DOG 
is less luminous than S 850 > 5mJy SMGs, some frac- 
tion of the DOGs are related to SMGs as shown in the 
30% of SMGs which meet the DOGs criteria. In order to 
obtain a submm sample of comparable number density 
to the 5*24 > 100 /iJy DOGs sample requires a survey 
down to 5s5o > 3mJy (Coppin et al. 2006). This will 
be achieved with future deep SCUBA-2 surveys and al- 
low for a more detailed comparison between DOGs and 
submm-cmitting galaxies. 

5. SUMMARY 

From a sample of 79 faint (S24 > 100 /iJy) DOGs in 
GOODS-N (0.5 arcmin -2 ), we find that almost all satisfy 
the criteria for Compton-thick AGN from F08. However, 
based on Spitzer spectroscopy and photometry, we show 
that 80% of are likely dominated by star formation. The 
stacked X-ray emission from the mid-IR classified star 
forming DOGs is consistent with what is expected from 
star formation. 

The IRS spectra and Spitzer photometric redshifts con- 
firm that these faint DOGs lie in a tight redshift distri- 
bution around z ~ 2. Stacking the mid-IR, far-IR and 
submm flux of the star forming DOGs, we derive an av- 
erage SED with Lir ~ lx 10 12 L Q , 8 times less luminous 
than most bright (Ssso > 5 mJy) SMGs. The composite 
SED of DOGs has a similar shape to that of SMGs in 
the far-IR (dust temperature of around 30 K) but has a 
higher mid-IR to far-IR flux ratio (Lir/Ls ~ 7 compared 
to Lyr/L% ~ 20 for SMGs). This suggests that there is 
a wide range of Lir/Ls conversions in z — 2 galaxies 
which need to be considered when interpreting the total 
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IR luminosity density and SFRD from 24 /im surveys. 

The average star forming DOG has a star formation 
rate of 200M Q yr~ 1 which amounts to a contribution 
of 0.01 Mq yr _1 Mpc" 3 (or 5-10%) to the star formation 
rate density at z ~ 2. 

This paper has relied strongly on stacking analysis to 
obtain average properties of DOGs (Fig.[SJ). Future deep 
surveys with Herschel and SCUBA-2 will detect the ma- 
jority of DOGs putting constraints on their individual IR 
luminosities and dust temperatures. 
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